Superconductivity in the Hubbard model with correlated hopping: Slave-boson study 
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The slave boson mean-field studies of the ground state of 
the Hubbard model with correlated hopping were performed. 
The approach qualitatively recovers the exact results for the 
case of the hopping integral t equal to the correlated hopping 
integral X. The phase diagram for the strongly correlated 
state with only singly occupied sites, the weakly correlated 
state, where single and double occupation is allowed, and for 
the superconducting state, was determined for any values of 
X and any electron concentration n. At the half-filled band 
(n — 1) a direct transition from the superconductor to the 
Mott insulator was found. In the region of strong correla- 
tions the superconducting solution is stable for n close to 1, 
in contrast to the case of weak correlations, in which super- 
conductivity occurs at n ~ and n ~ 2. We found also 
that strong correlations change characteristics of the super- 
conducting phase, e.g. the gap in the excitation spectrum has 
a nonexponential dependence close to the point of the phase 
transition. 

74.20.-z, 71.27. +a, 71.30.+h 



I. INTRODUCTION 

Over 10 years passed from the discovery in cuprates 
superconductivity with a high critical temperature, how- 
ever, its nature and the mechanism leading to super- 
conductivity have not been yet identified in a conclu- 
sive manner. There are many experimental facts, which 
suggest a significant role of magnetic and charge-charge 
interactions in formation of the superconducting phase 
(e.g. in cuprates and bismuthates, respectively). The 
simplest model describing the physical situation in such 
compounds and taking into account electronic correla- 
tions on the lattice is the Hubbard model. An important 
question, which arises is a mechanism of formation of 
Cooper pair in the presence of strong Coulomb onsite 
repulsion. In some circumstances even not strong inter- 
site Coulomb interactions, as that one for correlated hop- 
ping, overwhelm onsite reaulsion and lead to supercon- 
ductivity for hole hoping. pa\ Recently the ID Hubbard 
model with correlated hopping was analyzed for the case 
when the hopping integral t is equal to the correlated 
hopping integral X. [LTEJ] Exact results were obtained 
by using the conservation of the number of doubly oc- 
cupied sites in the model and mapping it onto a system 
of spinless free fermions with twisted boundary condi- 
tions [t|, or in a equivalent way, by expression of the 
model by slave bosons, construction of the eigenstates 



sy with respect of number of 
H| The superconducting phase 



and minimizing the en 
doubly occupied sites, 
was found as a ground state for the onsite Coulomb in- 
tegral U < —Atcos(nn) (n is the concentration of elec- 
trons). In the case of strong attraction, U < — 4t, in the 
ground state the sites are occupied by two electrons or 
are empty, and there are no singly occupied sites. In the 
region U > —At cos(7rn) there are only singly occupied 
sites allowed, what corresponds the situation U = oo in 
the pure Hubbard model (X = 0). 

The other type of intersite interactions may destabilize 
superconductivity. The density-density coulomb interac- 
tions lead to the charge density wave (CDW) ordering for 
n = 1. The stability conditions of the CDW were stud- 
ied in .an exact manner for.-n = 1 and X — t by Strack 
et al. pi] and Ovchinnikov [0]. The stability of saturated 
ferromagnetism in the presence of the exchange Heisen- 
berg and. the pair hopping terms was studied by Kollar 
et al. pj] using the generalized Nagaoka's theorem for 
n = 1 with one hole. The most general|-situation was 
analyzed by de Boer and Schadschneider f5\. Using the 
optimal ground state approach they determined [ES] suf- 
ficient conditions for the uniform and modulated super- 
conducting state, CDW, AF and ferromagnetic state for 
the extended Hubbard model, which included all type 
of nearest neighbor (n.n.) interactions on the arbitrary 
lattice at the case X = t. 

Deviation from the special point X = t removes the de- 
generacy of the ground state. There are not known exact 
analytical methods, which work apart from this point for 
arbitrary values of X, and may be use to determine the 
ground state. The Hartree-Fock approximation (HFA) 
may be useless for such studies as it neglects correlations 
of electrons and can not describe the Mott transition at 
n = 1 as well as the state with singly occupied sites (the 
strongly correlated state) for n ^ 1. The slave boson 
mean-field approximation pa] (SBMFA) may be use for 
this purpose as it is equivalent to the GutzwiHer approx- 
imation to the Gutzwiller wave function. p3] The ap- 
proach has been extensively usedr&s-the Hubbard model 
and its various generalizations. [tZrEZl] For the repulsive 
case U > the stability of different magnetic orderings: 
ferromagnetic, ferrimagnetic, antifcrromagnetic and spi- 
ral structures, as well as the CDW state were investi- 
gated. [EirEa] An influence of a static distortion of the 
lattice on— the stability of the magnetic phase was also 
studied, pa] Properties of the superconducting state were 
studied either by a transformation of the attractive Hub- 
bard model to its repulsive analog and solving the stabil- 
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ity conditions for the antiferromagnatic phase in the pres- 
ence of an external magnetic field |E9] , or by means of the 
spin- and the chaxgaj-otationally invariant repcesentation 
of slave bosons p§E3]. The SBMFA results p3] for the 
ground state (the free energy, the gap in the excitation 
spectrum, the critical magnetic field and the coherence 
length) were compared with those obtained in the HFA 
and by the self-consistent second order perturbation the- 
ory in the weak coupling regime (U/t <C 1) |Oc3] as well 
as with those obtained using perturbations approaches in 
the strong coupling regimes (U/rt-p> 1) [El]. The SBMFA, 
in contrast to the HFA, gives pJ] credible results in the 
whole interaction range and for any electron concentra- 
tion. The one orbital Hubbard model with a generalized 
correlated hopping term has been derived as an effective 
model from the three orbital model for Cu02 plane of 
cuprates by Simon and Aligia pH]. They also used the 
SBMFA, but confined themselves to the analyze of the 
Mott transition for some restricted range of parameters 
of the model. 

The purpose of the present work is to determine the 
ground state of the Hubbard model with correlated hop- 
ping and to study its basic characteristics in for whole 
range of parameters X and U, as well as for any elec- 
tron concentration n. We are convinced that the SBMFA 
gives reliable qualitative description of the physical sit- 
uation, especially the strongly correlated normal phase 
(with only singly occupied sites) and the superconduct- 
ing phase in the presence of correlations. After the intro- 
duction of the model and the method in the Section 2, we 
compare, in the Section 3a, the SBMFA results for X — t 
with those obtained in the HFA as well as with rigorous 
solutions. It will be shown that the SBMFA results are 
very close to the exact ones, and differ from the HFA. 
The case of the half-filled electronic band (n = 1) is con- 
sidered separately in the Section 3b. We will determine 
the phase diagram for the superconducting state in the 
space of the parameters of the model. Apart from the 
phase transition from the superconducting to the metal- 
lic state for small values of A, we expect to find the 
phase transition from the superconducting state directly 
to the Mott insulator in a finite range of X close to the 
point X = t. We will undertake investigations of char- 
acteristics of the superconductor-to-Mott insulator tran- 
sition, as due to strong correlation effects they may be 
different from those known for the superconductor-metal 
transition. The Section 3c is devoted to the analyze of 
the stability of the solutions for the superconducting and 
the normal phase in the most general case n ^ 1 and 
I^i. We predict the strongly correlated metallic phase 
around the point X = t. The region of the stability for 
this phase depends on the concentration of electrons n, 
and it can be different for doping by holes and electrons. 
The phase transition to the superconducting state will 
be studied in the whole range of X and n. We will com- 
pare the SBMFA results with those obtained by the HFA. 
The differences may be quite large, especially close to the 
point X — t, where the HFA does not work properly and 



can not reproduce the exact results. pilS] Discussion and 
comparison with other approaches will be performed the 
last Section. 



II. SLAVE-BOSON REPRESENTATION OF THE 
HUBBARD MODEL WITH CORRELATED 
HOPPING 

The Hamiltonian of the Hubbard model with corre- 
lated hopping is given by 

H = ( c i<T C j° +4 <r c ia )[-t + X(n i - <r +n j - .)] 

<i,j>,<J 

+U2_,n^n iU (1) 

i 

where the sum is restricted to n.n. sites. The param- 
eters of, the Hamiltonian (|l|) were determined by Hub- 
bard p3] as U — 20 eV and X — 0.5 eV from purely 
Coulombic origin. Screening effects may considerably re- 
duced their values, p3] especially in highly anizotropic 
systems. [E3] An estimation for cuprates and bismuthates 
gives U = 3 5 eV, X = 0.1 -h 0.3 eV and t = 0.1 ~ 0.2 
eV. There are several mechanisms, which may lead to a 
further reduction of U and even to overcome it. One of 
them is a polaronic mechanism, which can be realized in a 
narrow-band electronic system coupled to high-frequency 
intramolecular vibrations or cation-ligand vibrations. [El] 
The short-range attraction can result from a coupling 
between electrons and-cmasi-bosonic excitations such as 
excitons or plasmons. [Ej] Another possibility is a purely 
electronic origin resulting from a coupling of electrons 
in a many opbital model as has been derived by Simon 
and Aligia. p3\. Therefore, the Hamiltonian (P may be 
treated as an effective one, where its parameters U and 
X may change in the whole range, they may be also neg- 
ative. The hopping integral t is taken positive (without 
losing generality of the considerations). 

We want to express the Hamiltonian (|l|) in the slave- 
boson representation, Wrhieh is the spin- and the charge- 
rotationally invariant. [EaH] The singly occupied states 
| a) i are expressed by the Bose operators p\ , and the 
Fermi operators fj a as 

W)i = Y,A™>fU vac ) ( ff =t,i). (2) 

a' 

The doubly occupied state |2)j and the state |0)j, corre- 
sponding to an empty site, form a doublet (p — 2,0), 
which may be expressed by the Bose operators b\ , as 

p' 

where 
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ft f t 

1 



(4) 



In order to operate in the physical part of the extended 
Hilbert space we introduce the following constraint: 



2 E CAp'p + 2 E^U'P-'- = L 



(5) 



It ensures that each site is occupied by exactly one slave- 
boson. The Hubbard operators are given by 



\<J\)i i(0" 2 



2 E?U 



Pi. <j<j2 



x 7 = w)i *{p\ = E^l^'4^' E hp'p^p' 

a 1 p' 
= ^2pta*>(K0pft<T> - o'k,2pfi-v')< (6) 

By means of these operators any extended Hubbard 
Hamiltonian may be expressed, especially the Hamilto- 
nian (1), which takes the form 



i?= E [hXfX° a +t 2 X 2a X 



cr2 



<i,j>,a 



-ot l2 {x1- a xf + XfXJ* 2 ) 



h.c. 



(7) 



where £j = — t, f 2 = —t + 2X, t\ 2 = —t+X. For arbitrary 
values of the parameters t±, t 2 and t\2 it is the most 
general form of the Hamiltonian with correlated hopping. 
If t\2 = the Hamiltonian commutes with the total spin, 
the number of doubly occupied sites D = J2i n i1 n il — 
J2i Xf 2 and the r] operators: 



r? = E c ^i| = -E X ° 2 ' 

i % 



HI 



1) 



(8) 



This fact was y^ed in finding the ground state in the 
exact manner. pi3] 

We want to study the Hamiltonian (]?]) in the mean field 
approximation, but instead simply to exchange the Bose 
operators in Eqs. (|6]) by the c- numbers we renormalize 
properly the operators X. We use the fact that repre- 
sentation of any operator by the slave operators is not 
uniquely defined, but only up to operator factor whose 
eigenvalue is unity in the physical space. The renormal- 
ized operators X may be expressed as 



X°» = 4»fl-aA¥n^, 



(9) 



where the 2x2 matrix A is 

A = v'LRb, 

and 

L= [l-pt £ -6tg|-V2 , R= [l-ptp-^b]- 1 / 2 
b 



(10) 



JQ2 "00 



b 22 b 2 Q \ g _ / On,, - - l) 2 n 



boo -b 2 
—bo2 b 22 



„ = . PTT PU \ g = ( PU -PTi 

\pit pu/ ' - V^PiT PTT 



(11) 



Now, we may perform the mean-field studies of the 
stability of the superconducting phase. We determine 
the lowest energy solution from the minimum of the free 
energy F 

F = F f + F b = 

~\ E^ 1 + cx P(^)} + ln il + exp(-^ fe )}] 

+ — (& 2 + 2,5)-Ao(l + 2<5)iV-2AAAiV, (12) 

where Ff and -Ff, are the fermionic and bosonic parts of 
the free energy, A is the superconducting order parame- 
ter defined by 



A = (c uCiT > =2^(&t p2 6 i0p ), 



(13) 



and E k = sj {-2tr\k + Xo) 2 + (2rrjk + Aa) 2 is the disper- 
sion curve of fermions described by the Hamiltonian (0) 
in the SBMFA, i.e. by 



k a 

+ (2r Vk + \ A )(fl ] fl l +h.c. 
The following notation is used 



(14) 



I = -a 2 ^! +t 2 + 2*i 2 ) + 2Sb 2 {-t 1 + t 2 ) 

+ {5 2 - A 2 )^ + t 2 - 2t 12 )}, 
r = a 2 [b 2 (-h + t 2 ) + S{h +t 2 - 2t 12 )]2A, (15) 

where 

a 2 = P 2 /[6 2 (l-4£ 2 -4A 2 )], 



6 2 = (6 2 + v 7 ^ 4 -4<5 2 -4A)/2 , 

&2 = 2 < & 1 2 2^22 + &lo(AoO + &J 20 & l0 2 + &| 2 & i20) , 

P 2 = 2 <PiTTPiTT +pl U Piil) , 
the electron concentration is 

n = 1 + 25 = 1 + 2((4 2 6 l22 ) - (0*00)) , (16) 

r/fc = J2 a cos ^a, -W is the number of the lattice sites. The 
Lagrange multipliers Ao and Aa are for the constraints, 
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which ensure the number of electrons and the supercon- 
ducting order parameter A. 

The free energy (|l2|) is the function of the variational 
parameters: 6 2 , A, Ao and Aa- The parameter b 2 changes 
in the range from 2\ / S 2 + A 2 to 1, and the s upercond uct- 
ing order parameter A < i yl — 4S 2 — | y/n{2 — n). In 
the normal phase (A = 0) b 2 describes the concentra- 
tion of doubly occupied sites, which is if fo 2 = 2\S\, 
and maximal equal to n/2 for b 2 = 1 (according to nota- 
tion used by Kotliar and Ruckenstein |t£j in the normal 



phase b 



2d + 1 — n, where d and e are the 



concentration of doubly occupied sites and empty sites, 
respectively). At b 2 — 2\5\ there are no doubly occupied 
sites, the state is paramagnetic metal with a maximal 
value of the local magnetic moment. Such a state cor- 
responds the situation for U = oo in the pure Hubbard 
model (the case X = 0) and we call it the strongly cor- 
related normal state (SC). For 2|<5| < b 2 < 1 there are 
allowed doubly as well as singly occupied sites. We call 
it the weakly correlated normal state (WC) in order to 
differentiate it from the SC state. The local magnetic 
moment decreases with increasing b 2 , and disappears at 
all for b 2 = 1, i.e. when only doubly occupied sites are 
allowed. For ther-nure Hubbard model (A=0) it happens 
in the SBMFA |E3] for strong attraction (U < 0), when 
the effective bandwidth becomes zero, the system is then 
a diamagnetic insulator. Comparison of the SBMFA free 
energy with that one obtained by the perturbation ex- 
pansion in the limit |Z7| 3> t shows p3] that there is still 
a part corresponding to mobile particles in a band, what 
neglects the SBMFA. For large \U\/t the Cooper pairs 
have very small size, which is restricted to a single lat- 
tice site, and they form a Bose liquid. The system in the 
normal phase is a diamagnetic metal, in which electronic 
transport is due to motion of these local pairs (hard core 
bosons with charge 2e). Condensations of the bosons 
leads to superconductivity [El]. 



III. THE SBMFA GROUND STATE 
CHARACTERISTICS 

A. The t = X Case 

Here, we want to analyze the ground state solutions 
in the SBMFA. In order to see validity of the present 
approach we perform calculations for X = t and compare 
them with -the results obtained in the exact way for the 
ID chain p]. The considered case will be a reference 
system for further studies. It is worth to notice that 
in high temperature superconductors the parameters are 
estimated as X = 0.1-f-0.3 eV and t = 0.1-^0.2 eV, what 
is close the considered case. 

At X — t one has an electron-hole symmetry, which in 
general (for t^ X) is not kept. The parameters of the 
Hamiltonian (|l4|) are now 



t = -4a 2 6gi(5 and r = Aa 2 bltA. 



(17) 



To determine the free energy ( |l2| ) we first transform the 
variables to 



S = (tS + rA)/p, A = (r5 - tA)/p, 
A = (t\ + rXA)/(2p 2 ) , X s = (A r - 



X A t)/(2p 2 ), (18) 



where p 2 = 



It corresponds to a rotation in the 



space of the parameters S and A. The dispersion curve 

is expressed as Ek = 2p\J (ijk + Ao) 2 + Ag. For t and r 
given by (0) the parameter A = 0, what implies A s = 
as well. The gap in the excitation spectrum, defined as 
Eg = 2max(£^), is also equal to zero in this case. The 
free energy of the superconducting state for the ID chain 
at T = is 



F s /N = -2p- 



dk\r, k + X Q \ + j(b 2 + 2S)-4\ Sp. 

(19) 



The minimum of F is at Ao = —25, b 2 — 2\S\ and the 
solution may be expressed in the parametric way as a 
function of 6 and \S\ = VS 2 + A 2 (\S\ < \5\ < 1/2) 



F s /N 



16 ^l .sm^(l-2\8\) + U(\S\+5), 



tt(1 + 2|£|) 2 



(20) 



c/= ^sinf(l-y|) , 
1 + 2|<5| L 7r(l + 2|S|) 11 2 V 



(21) 



In a similar way we determine the ground state energy 
for the normal phase (A = 0). There are only two stable 
solutions: with b 2 = 2\5\, corresponding to the state of 
strongly correlated electrons, in which only single occu- 
pancy is allowed, and with b 2 — 1, corresponding to the 
diamagnetic metal, in which single occupancy is forbid- 
den. There is not present the weakly correlated phase 
(with 2 1 5 1 < b 2 < 1), in which occurs sites with two elec- 
trons as well as sites with one electron. The free energies 
for these two stable situations are 



Ff/N = - 



F»/N = U~ 



for U>U 12 , (22) 
for U <U 12 , (23) 



where 



U 



12 



m\s\ 



sin — (1 



7r(l-45 2 ) 2 



2\S\) 



(24) 



Putting A = (\S\ = \S\) in Eq. @ we get the criti- 
cal value of U, below which the superconducting phase 
is stable. The phase transition is the second order to the 
normal phase of strongly correlated electrons. The ex- 
citation spectrum in the superconducting phase has no 
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gap, however, there are stable Cooper pairs (A ^ 0) as it 
is seen from the energy difference F N — F shown in Fig. 
|l|. For U < — At the free energy of the superconducting 
state and the normal state are the same. This degener- 
ation may be removed by any interaction stabilizing one 
of the phases. We will show later that a deviation from 
X = t stabilizes superconductivity. 

The solutions in the Hartree-Fock approximation may 
be also presented in the parametric form with the param- 
eter s. The superconducting state has the free energy at 
T = 

Fhf/N = sin |(1 - 2s) + Us 2 + ^(1 +45), (25) 

U = At\— sin^(l - 2s) - cos^(l - 2s)} . (26) 

7TS 2 2 

and the order parameters are given by 



A = < (kicii >= y ■ 
Ai = < c a c jT >= 



5 2 



s 

-sm-(l-2 S ) 

7T Z 



(27) 
(28) 



This state is stable for U lower than 



U. 



HF 



4 ^ Sin f (1 ~ 2|,5|) " COS l (1 ~ 2|,5|)] - (29) 



The free energy of the normal state is that given by 
Eq.© with s = \S\. 

For comparison the exact results [&El] for the ground 
state energy in the ID chain are 

F ex /N = --L v /l6f2_[ / 2 + ^ [n _ I arccos( _^)] 

for - At < U < Ul x , (30) 



2t 

F ex /N = sin 7m for U > Uf , 

7T 



Fex/N = U- 



for U < -At , 



(31) 
(32) 



where UJ? = — At cos im. 

Fig. g shows the quantitative comparison of the free 
energies determined by these methods for n = 1. The 
SBMFA solution (the solid curve) is very close to the ex- 
act one (the dashed curve) . Although the SBMFA and 
the HFA solutions [Eqs.@-@ and Eqs.©-©, re- 
spectively] are similar in their structure, the values of the 
free energies are different. The free energy for the HFA 
(the long-short curve in Fig. ||) is always higher than the 
exact and the SBMFA ones. The reason is in a reduction 
of the hopping integral and the width of the electronic 
band (the narrowing band effect) in the SBMFA. The 
physics in the HFA is also different: in the normal phase 
may occur doubly occupied sites even for very large U. 
For U < —At the system is highly degenerate, and all so- 
lutions [Eqs.©, @, © and '(H)] (as well as the free 
energy for CDW at n = 1) are identical F/N = Un/2. It 
is energy of localized electron pairs, in an absence of hop- 
ping process, and therefore, it corresponds-ito an insulat- 
ing state. The Drude weight determined pi] by changing 



the twisted boundary conditions of the ID model equiv- 
alent to (0) (and similar to (0)) is D c = for U < -At 
and any n, indicating on the insulating phase. Above 
Us the stable solution is the strongly correlated state. 
The same free energy has the ferromagnetic state with 
saturated magnetic moment m = n/2 (as well as the 
antiferromagnetic state in the case n = 1). In the in- 
termediate regime —At < U < Us, may also occur a 
nonuniform ordering, in particular, the state with two 
separated phases: the phase with only doubly occupied 
sites and the phase with only singly occupied sites. The 
stability of a such two-domain system may be analyzed 
by the Maxwell construction, using the single domain re- 
sults for the SBMFA, Eqs.@ and @. The free energy 
is then exactly the same as for the uniform superconduct- 
ing state given by Eq.(p0|). Onegan show that the results 
obtained in the exact manner p a [E qs.(|30|)-(|32])1 have 
the same properties (see also Ref . [ B6[ ) . 

The SBMFA may be apply also for the systems of 
higher dimensions. The results depend on the density 
of states (DOS) and differ slightly. For illustration we 
present below the ground state characteristics for the 
rectangular DOS (N(E) = l/W for \E\ < W/2): 

F s /N = -^(1- ^- f + US for -W<U<U S , 



F s /N = U- for U < -W , 
F»/N = -W\S\{1 - 2\S\) + U(\5\ + S) 



for U > -2\S\W , 



FJ?/N = U- for U < -2\6\W 



(33) 



(34) 



where Us — W(1 — A\S\). These results may be compared 
with those analyzed above as well as with those for a 
d-dimensional hypercubic lattice taking the bandwidth 
W = Adt. The HFA applied to the Hamiltonian (0) with 
the rectangular DOS gives qualitatively different results: 



U. 
A 



F* F /N=^(l + A6) 



1 



108VF 2 



[U - VU 2 +3W 2 ][6W 2 + U(U - VU 2 + 3W 2 )} 



for — W < U < U§ F , 



U 



Fhf/N = -n for U < -W , 
F H T F /N= ^n 2 -^\S\(l-AS 2 ), 



(35) 
(36) 

(37) 



and C/f F = W(l -12|£| 2 )/(4|(5|). 

Fig. U summarizes our results. There are exhibited 
the dependencies of Us vs. the electron concentration n 
for the different approaches. The HFA curves (the long 
dashed curve and the curve with crosses) are completely 
different than the others. Of course, in the HFA there 
is no place for the strongly correlated normal state with 
only singly occupied sites. With an increasing deviation 
from the point n = 1 (increasing \n — 1|) all the curves in 
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Fig. H become closer to each other, correlations between 
electrons become weaker. 



B. The Mott Transition in the Hubbard Model with 
Correlated Hopping: The n = 1 Case 

The case of the half- filled band (n = 1) we would like 
to treat separately as it is a discontinuity point of the 
model. Moreover, for X = (in the Hubbard model 
without correlated hopping) there is the Mott transition 
as a function of U from a metallic phase to an insulat- 
ing phase, in which double site occupancy is forbidden. 
We have already learnt (in the previous section) that at 
X = t the strongly correlated state may occurs for any n 
and it is metallic apart from the case n — 1 correspond- 
ing to an insulator. The Mott transition is then from the 
superconductor to the insulator. The exact diagonaliza- 
tion calculations for the ID chain performed by Arrachea 
et al. |E3] at X = — t indicate also such a transition. 

Our procedure is similar as in the previous section. We 
perform our further calculations for the rectangular DOS, 
which are much simpler and give qualitatively similar 
results as for DOS of d-dimensional hypercubic lattices. 
After the transformation ( |l8| ) the free energy at T = is 



b 2 

Uj-4p(\ 6- 



A S A) 



where 



f s = d{R + + R_) + \ {R+ - R-) 

-2 R + + X Q + d 
+Ao In = , 

6 R- + X Q - d 



R± = V(A ±d) 2 + Ac 



(39) 
(40) 



The bandwidth of noninteracting electrons is W = Adt. 

First, we analyze the normal phase (A = 0). There 
are possible three type of solutions, in which: (i) double 
occupancy is forbidden (b 2 — 0) for large U - the Mott in- 
sulator, (ii) there are sites with two electrons as well sites 
with one electron (and, of course, there are also empty 
sites) - in an intermediate coupling regime , (iii) single 
occupancy is forbidden (b 2 = 1) - in a strong attractive 
regime. We find that the minimum of the free energy 
(E38|)-(E(t) in the normal phase is: 



F N /N = 
F N /N = -d\t - X\ 



for U > Umo 



U 



1? 



F N /N = 



U 



8d\t-X\ 

for \U\ < Umou , 

for U < -Umou , 



(41) 

(42) 
(43) 



where the critical value for the Mott transition is Umou — 
8d\t — X\. It is seen that the results are the same as in 
the absence of the correlated hopping (X = 0) if we use a 



renormalized hopping integral \t — X\. The average con- 
centration of doubly occupied sites is b 2 = |(1 — ggrjzyi ) 
and monotonicly increases from b 2 — at U = Uuott and 
reaches b 2 — 1 at U = — Umou- 

The analyze of the superconducting phase is now more 
complicated as we have to find the minimum of F, 
Eq.(p8[)-(^0|), which is a function of the parameters b 2 , 
A for a given U and n = 1 + 26. If we assume that the 
transition to the superconducting state is the second type 
we can expand F in a series of A and find Us - the critical 
value of the parameter U below which superconductivity 
is stable . There are two solutions corresponding to the 
transition from the metallic phase: 



U s = 8d(\t -X\- x/t 2 - 2tX + X 2 /2) 



and from the insulator 



U s = 4dy/t 2 - 2tX + 2X 2 



(44) 



(45) 



The stability diagram for the superconducting phase is 
presented in Fig. ||, where the solution Q) and fl45| ) 
are shown by the solid curves. The dashed line shows 
the stability regions in the normal phase. In the region 
of small X the phase transition to the superconducting 
state is from the metallic state; for larger values of X the 
transition is direct from the Mott insulator. For compar- 
( 38 ) ison, the critical value U¥ F = 2dX 2 /\t - X\ obtained in 
the HFA is shown in Fig.[| by the long-short dashed curve 
(see Appendix C for details). The significant differences 
between the both approaches are for X > 0.5t. 

Our studies do not confirm the prediction by Arrachea 
et al. p3] the superconductor-to-Mott insulator transi- 
tion at X = —t. They performed exact diagonalization 
calculations for the ID chain, in which the Mott insula- 
tor is the stable solution for any positive U. An compar- 
ison with exact results obtained [til] by the Bethe Ansatz 
for the attractive Hubbard model (at X = 0) shows that 
some fluctuations specific for the ID system are not taken 
into account in the SBMFA. p3\ The superconductor-to- 
Mott insulator transition in ihe ID chain for X < 0, 
obtained by Arrachea et al. |E3], may not occur in sys- 
tems of higher dimension, where one expects the Mott 
transition for a finite value of U and the SBMFA works 
better. 

The stable solutions and their ground state character- 
istics are, in general, determined numerically. In FigJEja 
we present the free energy difference AF = F N — F s 
between the normal and the superconducting state vs. 
the parameter U for a few given values of X. The curves 
for X — 0.8t, 1.2t and 1.5t correspond to the solutions 
of the superconducting state in the strongly correlated 
regime, where double occupancy is forbidden. The AF 
has a maximum at X = for all these cases. We see that 
at X = t the free energy for the superconducting and the 
normal phase is the same for large attraction (U < —W), 
but any deviation from X = t favors the superconduct- 
ing solution (AF > for any U < Us). At X = 0.5* 



G 



and X = — 0.8t the system is in the weakly correlated 
regime, where the metal-to-superconductor transition oc- 
curs. The maximum of AF is now at ca. 2|1 — X/t\W 
(the dotted curve has its maximum outside the figure at 
U = — 3.1W), similarly as in the pure Hubbard model 
(X = 0) (see Ref.[ for details). The U dependence of 
the gap in the excitation spectrum, E g = 4p\Xs\, is also 
different in the both regions, ft is visualized in Fig.||b. 
For the strongly correlated case (see Appendix A for an 
analytical solution) the effective bandwidth W e tf = Adp 
as well as the gap in the excitation spectrum E g linearly 
increases with Us — U (see Fig.||c and [5(d). Therefore, the 
value of Eg is relative small in a wide region of U. In the 
weakly correlated region and for U close to Us the en- 
ergy gap Eg has an exponential dependence. It is similar 
to the Hartree-Fock solution, but the value is reduced by 
the factor 7 = 0.47. |E3] 

In Fig||a the difference of the free energies of the nor- 
mal and the superconducting phase is shown as a func- 
tion of X for different values of U. The relative sta- 
bility of the superconducting phase is a nonmonotonic 
function and its character may be quite different for the 
different values of U. The dependence of the gap E g in 
the excitation spectrum is shown in Figj^b. E g has the 
maximum in the weakly correlated region, reaches zero 
at X = t and monotonicly increases for larger X. The 
small jump around X w 0.65i is at the crossover point 
from the region of weakly correlated electrons to the re- 
gion of strongly correlated electrons. The dependence of 
W e ff on X is shown in FigJ|c. 

Magnetic properties of the Hubbard model with cor- 
related hopping are apart from our main interest in the 
present paper, but at n = 1 the magnetic ordering may 
significantly influence superconductivity, especially for a 
positive value of U. Here, we want to present our stud- 
ies of the relative stability of the antiferromagnetic (AF) 
and the superconducting state in the considered model 
for n — 1. To analyze the stability of the AF and the 
superconducting state is needed to compare their free 
energies at a given point of the parameters of the model. 
The procedure of calculations of the free energy of the 
AF state is sketched in the Appendix D. It is simpler 
to determine than for the superconducting case, because 
one determines the free energy only for X = and the 
proper rescaling of the bandwidth gives the values for 
any X. Fig.]?] shows the phase diagram. At X = t the 
transition between both phases is at U/W = 1. The local 
magnetization is then maximal m = 1/2. For small X 
(X < 0.66*) as well as for large X (X > 2.08t) the bound- 
ary line separates the AF phase from weakly correlated 
superconducting state, whereas in the intermediate range 
- from the strongly correlated superconducting state. We 
also analyzed the relative stability of the both phases in 
the HFA. Their free energies are always higher than those 
obtained by the SBMFA. The HFA boundary is shown 
by the long-short dashed curve in FigJ^. It is seen that 
for X > it lies much lower than the SBMFA line. One 
can say that correlations (included in the SBMFA) more 



prefer superconductivity than magnetic ordering. The 
CDW phase does not occur in FigJ^ as its free energy is 
higher than that one for the AF and the superconducting 
state in the whole range of the parameters (apart from 
the case X = t and U < —W as well as X = and U < 0, 
where the ground state is degenerate). 



C. The case of arbitrary X and n 

Let us now analyze the ground state properties of the 
model apart from the special points {X — t and n — 
1) discussed above. Although the present case is more 
complicated, we may follow in the way described already 
in the Section 3a and 3b. The effective hopping integrals 
of the Hamiltonian (|l4|) are 



(1 - b 2 )[(X - t){b 2 + V6 4 - 4A 2 - AS 2 ) + 2X8} 



1 - 4A 2 - AS 2 



l-b 2 



1 - 4A 2 - AS 2 



2AX 



(46) 



It is also needed to transform the variables S, A, Ao and 
As according with the formula ([l8|). We perform the 
calculations for the rectangular DOS, for which the free 
energy at T = is expressed as 

e p Q b 2 + 25 4p ,~ - ~ — 

F /N = £-/ s + U 



where 



2d J 



-^(X 5 + X S A), (47) 



T 2 , 

+Aeln 



R-) + \ (R+ - R_) 
R + + X Q + d 



R-+X 



R± = \^(Xo±d) 2 + X s 



(48) 
(49) 



We study first the stability of the normal phase (A = 
0). Our task is to find, for given values of the parameters 
X, U and the electron concentration n = 1 + 28, the 
minimum of the free energy 



F N /N = -2d(l - b 2 )\{X - t)(b 2 + y/b* - AS 2 ) + 2XS\ 
+ j(b 2 + 25), (50) 

where the variational parameters b 2 is from the interval 
[2\S\, 1]. By comparison of the free energies we find the 
stability diagram for the normal metallic state. The sit- 
uation for electrons (n < 1) and for holes (n > 1) is 
presented in Fig.||a and Fig.^b, respectively. The sta- 
ble solution with b 2 — 2\8\ corresponds to the state of 
strongly correlated (SC) electrons (holes) and it occurs 
(above the solid curve) for X > t when n < 1 or for 
0.5t < X < t when n > 1. The lowest solid curves in the 
figures correspond to the solutions obtained for n — > 
and n —> 2, respectively. It is seen that the boundary 
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line for n = 0.99 and n = 1.01 (the upper curves) are 
different than that one for n = 1 shown by the dashed 
line in Fig.||. The case n = 1 is a discontinuity point, and 
has been already considered. The stable solutions for the 
diamagnetic state with b 2 = 1 (DM) are in the region 
below the dotted curve. The weakly correlated metal is 
stable in the region between the solid and dotted curves. 

The normal state is characterized by the parameter b 2 . 
We show its dependence on X in Fig.^a and Fig||b, for 
n = 0.5 and n = 1.5, respectively. There are two types 
of dependencies. If U is large enough an increase of X 
leads to a monotonic decrease of b 2 (see the lower curve 
in Fig.^Ja), which reaches the value 2|<5| at the boundary 
with the region of strongly correlated electrons. With 
a further increase of X the value of b 2 is constant up 
to the transition to the region of weakly correlated elec- 
trons, where b 2 > 2\S\. There is a first-order transition at 
the boundary line. As the upper curve in Fig.^Ja presents, 
the system may become diamagnetic (b 2 = 1), if the po- 
tential U is attractive and strong enough. One can show 
that the external magnetic field aligns all spins in the 
strongly correlated normal phase. The HFA calculations 
for the ferromagnetic state (Appendix C) give saturated 
ferromagnetism with the same value of the free energy as 
for the strongly correlated state. The HFA stability di- 
agram for the ferromagnetic state [Eq.( |C12] )1 is different 
than that one in Fig.|| as the free energies of the normal 
metallic phase determined by the both method are differ- 
ent. For holes (Fig.|^b) the situation is identical under the 
transformation t — > — t + 2X (the curves in Fig||b have 
the same character as those in Fig.^Ja if X decreases). 

Correlated hopping changes the width of electronic 
band and, in consequence, an effective mass of elec- 
trons. FigJ^c and||d show the effective bandwidth W e ff 
vs. X for electrons (n = 0.5) and holes (n = 1.5). 
W e ff is strongly dependent on X in the weakly corre- 
lated phase. Around the point X = t, W e ff is constant 
and equal to W/(l + \n— 1|) in the strongly correlated 
normal state or to in the diamagnetic state, respec- 
tively. The HFA effective bandwidth is simply given by 
WM = |1- Xn/t\W. 

Stability of the superconducting phase is determined 
numerically by comparison of the free energies for the 
state of strongly correlated electrons (b 2 — y/ S 2 + A 2 ) 
and weakly correlated electrons (the variational param- 
eter b 2 > \/S 2 + A 2 ). The phase transition from the 
superconducting to the normal state is a second-order. 
The stability diagram for the superconducting state is 
obtained from solutions in the limit A — > 0. It is pre- 
sented by solid curves in Fig.^Oja and [l0| b for different 
values of n < 1 and n > 1, . The dashed curves in 
Fig.[l(^ denote the regions of the normal phase of strongly 
correlated electrons (the same as in Fig.|J). The curves 
corresponding to n — and n — 2 are the exact solutions 
obtained for the case of two electrons and two holes in the 
Hubbard model with correlated hopping (see Appendix 
B for details). The SBMFA solutions converges to these 



ones if n — > and n — > 2. The point X — t/2 is singular 
for the solutions with holes. For the case of two holes the 
Cooper pair is stable for any value of U. In the SBMFA 
one gets at this point the superconducting phase stable 
for any finite concentration of holes (n > 1) and also for 
any U. 

Comparison, the critical values Us obtai ned by means 
of the SBMFA and the HFA [given by Eq.(^ is shown 
in Fig.|ll|a and |Tl| b for different concentration of electrons 
and holes, respectively. The SBMFA and the HFA results 
are closer to each other for n — > and n — > 2, as corre- 
lations become weaker and weaker. Small differences are 
also seen in the region X < O.bt. Strong correlations 
considerably modify the stability diagram for X > 0.5t, 
the results for the SBMFA and the HFA are qualitatively 
different there. 

The SBMFA studies for X = have shown |E3] that 
the value of the gap E g in the excitation spectrum can be 
much lower than the HFA result, especially for U close to 
Us- We expect a similar reduction of E g in the presence 
of correlated hopping in the weakly correlated region. 
Our SBMFA calculations confirm it. Fig.|l2| presents de- 
pendencies of Eg vs. n for X = — Q.8t (on the left hand 
side) and X = OAt (on the right hand side). The SBMFA 
and the HFA results have similar character, with a max- 
imum at small concentration of electrons (holes). 

The situation for X > 0.5< is different, as it shows 
Fig.jl^a. The superconducting state is stable for n close 
to 1, in contrast to the case X < O.bt from Fig.[l2[ There 
are the weakly correlated and the strongly correlated su- 
perconducting states for n > 1 and n < 1, respectively. 
The difference of the free energies of the normal and the 
superconducting state, F N — F s , exhibited in Fig.[l3|b 
indicates, that the stability of the both phases depends 
in a different way on X and U. For large U the stable 
solution may be only the strongly correlated supercon- 
ducting state. It is shown by the dashed curve, which 
corresponds to X — 3t and U — AW. Taking t = 0.1 eV 
we get X = 0.3 eV and U = 4.8 eV for a simple cubic lat- 
tice. Such the parameters are quite realistic for cuprates. 
The difference F N ~ F s as a measure of the condensation 
energy of the superconducting state is proportional to the 
critical temperature T c . For these parameters we esti- 
mate the value of T c of the order of 100K. From Fig.^ we 
see that at n — 1, X — 3t and U = WV the stable solution 
is the AF state, which should be also placed in Fig.[l3]. It 
resembles the phase diagram found for cuprates, where 
the AF state is destroyed and superconductivity emerges 
with a decreasing number of electrons in the system. 

IV. DISCUSSION AND FINAL REMARKS 

In this paper the role of the term of correlated hopping 
added to the Hubbard Hamiltonian has been studied by 
means of the SBMFA for the temperature T = 0. We 
have shown that the SBMFA stability of the supercon- 
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ducting phase is qualitatively similar to that determined 
in the exact way for X = t. This case is highly degen- 
erate as electronic correlations reduce the width of the 
band to zero for U < —W. We found the strongly corre- 
lated phase (with only singly occupied sites) as a stable 
solution also apart from the point X = t, namely, for 
n < 1 at X > t, and for n > 1 at t/2 < X < t if the on- 
site interaction U is strong enough. This state is metallic, 
in which all spins may be polarized in the presence of an 
infinitesimal external magnetic field. It means saturated 
ferromagnetism. The SBMFA stability boundary for this 
state is lower than the value of U determined by the op- 
timal ground state approach as thc-sufficient conditions 
for stable ferromagnetic solutions [ta]. At the half-filled 
band the strongly correlated state is the Mott insulator, 
however, it is not the state with the lowest energy. The 
ground state is either the AF state or the superconduct- 
ing state, depending on the value of U . The phase dia- 
gram for the superconducting phase is close to the HFA 
one only for negative and small positive values of X at 
any n, and for arbitrary values of X at n ~ or n ~ 2. 
For X < the superconducting state is stable at small 
concentration of electrons, whereas for < X < 0.5t 
at small concentration of holes. In the region X > 0.5t 
there are phases of the strongly correlated type. The 
stable solutions for the weakly and the strongly corre- 
lated superconducting phase occur close to n — 1. The 
U dependence of the gap in the excitation spectrum E g 
is also different there, it exhibits a nonexponential char- 
acter. For the parameters corresponding to cuprates the 
estimated value of T c is of the order of 100K. Our phase 
diagram, with the AF phase at n = 1 and the supercon- 
ducting phase for hole doping (n < 1), resembles that 
one obtained in cuprates. n 

The model analyzed by Hirsch and Marsiglio [□] differs 
slightly from that one considered above. They performed 
the particle-hole transformation of the Hamiltonian (|l|) 
with simultaneously done interchange k = and k = Q 
in the wavevector space k, i.e. c\ a — > (— l) 4 Ci CT , after 
which t —> t - 2X and X — > —X (see also Ref. [ || for 
the analyze of symmetry properties of the model) . In^ut 
ther considerations the hopping integral was taken 
as equal to t and independent on X . Their studies [ 
correspond to the case X < in the present work, where 
the weakly correlated superconducting phase occurs. The 
SBMFA boundary curve for the stability of the supercon- 
ducting state is a little bit higher than the HFA result 
(see Fig.[tl]). It is in agreement with exact dia«»aliza- 
tion calculations showing also a small increase. |E3e3] The 
SBMFA value of the energy gap E g is smaller than the 
HFA one. The same reduction is for the difference of the 
free energies, and consequently, for the thermodynami- 
cal critical magnetic field H c . In the strong attraction 
regime we expect superconductivity of local pairs (con- 
densed bosons). The coherence length £cl is then of the 
order of the lattice constant, and the H c decreases with 
increasing Us — U in contrast to the HFA, in which the 
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Hf? F always increases. .— . 

The Hubbard I approximationJu3] applied to the model 
was-nerformed by Das and Das [till] , and by Domahski et 
al. [H. In this approach two Hubbard s ub b an ds are sepa- 
rated by a gap for any value of U. They pd'cS] determined 
the critical temperature T c for the superconducting phase 
as a function of the parameters of the model and the 
concentration of carriers n. Although our analyze is re- 
stricted to T — we may compare the critical value of 
n, below which, at given U and X, the T c becomes posi- 
tive. Their resu lts a re in contrast to n c determined in the 
HFA [from Eq.(|C7|)] as well as in the present approach 
(the SBMFA) for the same parameters of the model. The 
reason of the discrepancy is in the method they used, 
which completely neglects interactions between electrons 
in the subbands even on the mean-field level. 

Our studies are restricted to the Hubbard model ex- 
tended only for correlated hopping, neglecting other n.n. 
terms. A most important one is the density-density 
coulomb interaction, which if it is repulsive may restricts 
the region of the strongly correlated normal state (see 
for example the exact analyze pEl] for X = t) and the 
superconducting state. The SBMFA correctly describes 
onsite electron correlations, but it is too poor for the 
n.n. density-density interactions. In the representation 
@"@ such the term is expressed by the operators b and 
in the mean-field approximation charge-charge correla- 
tions are properly taken into account (for example the 
CDW state [M]), whereas superconducting n.n. correla- 
tions are completely neglected. Some improvement of the 
present SBMFA is needed. We also confined our studies 
to the case T = 0, where fluctuations of bosonic fields 
were neglected. One can extend the studies for T > 
using the path integral approach and taking into account 
the fluctuations of bosonic fields in a coherent potential 
approximation (CPA) (in an analogous way as one treats 
static f3\ or dynamic fluctuations for D — oo pa] ). 
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APPENDIX A: SUPERCONDUCTING 
SOLUTIONS IN STRONGLY CORRELATED 
REGION FOR HALF-FILLED BAND AND 
CLOSE TO TRANSITION 

Here, we want to present analytical solutions for the 
superconducting state in strongly correlated regime for 
n = 1. In this case b 2 — 2|A|. Expanding the effective 
Hamiltonian parameters i, r [Eq.(|l5|)] in a series of A, 
we find that the effective hopping p — Apo, where po — 
Vt 2 - 2tX + 2X 2 , and from the condition dF/dA = 
the gap parameter A is found to be 



where 



A 



V-) 

4 l U s h 



(Al) 



where Us — AdpQ. The linear dependence of A is over 
a wide range of U for X close to t. It means that the 
effective bandwidth W e ff = Adp linearly increases with 
a decrease of U. From the conditions dF/d\o = and 
dF/dXs — we find in the limit As — > the expressions 
for U 



u = Us{1 + d[t-x-\ 0g Ts ) 

and for the gap in the excitation spectrum 



E 9 = M^\ = - ¥ ^ Y X s \o g ^ 



2 Po T 2 i„„^ 
As 



(A2) 



(A3) 



The other quantities may be expressed by means the 
above relations with As as the parameter. 



APPENDIX B: BINDING ENERGY OF TWO 
ELECTRONS (HOLES) IN THE HUBBARD 
MODEL WITH CORRELATED HOPPING 

The problem of two particles in an empty lattice may 
be solved exactly for our model. It is a good test for 
our calculations in the SBMFA. For two electrons on 
an empty lattice the trial function for formation of the 
Cooper pair is (see also Ref . [ fl3lj40| ) 



* = J2f kC h c -ki\ vac ) 



(Bl) 



From the Schrodinger equation for the Hamiltonian (|l|) 
and the trial function (Bl) we get 



(B2) 



The eigenvalues are determined from the condition 

1 - 8Xh - UI + 16X 2 (I 2 - I h) = , (B3) 



N ^ E + Atrjk 



(B4) 



We change the sum over k to the integral and calculate 
it under the assumption of the rectangular DOS (with 
\Vk\ ^ d) the corresponding integrals. The condition for 
the eigenvalue is now 

8dt(t X f = [Ut 2 EX(2t X)} log !^|| . (B5) 

The binding energy of the Cooper pair is Ef, = E — 4dt. 
From this condition we find the critical value of U , below 
which the Cooper pair is stable 



El 

8dt 



X(2t-X) 
21 2 



(B6) 



Changing t 
the system 



— t + 2X we get Us for the two holes in 
U s X(2t - 3X) 



8dt 2t\t-2X\ 



(B7) 



At X = t/2 the solution is singular, as the hopping for 
holes is equal zero. The binding energy of the Cooper 
pair is then 



E b = -(U - VU 2 + 6AX 2 a 2 ) 



(B8) 



where a 2 — jj J2k Vk ari ^ f° r the rectangular DOS 
a-2 = d 2 /3. For the ID chain the crit ical values of U 
are the same as those given by Eq.(p36|)-(B7) with d = 1, 



although the eigenvalues have different dependencies. 

APPENDIX C: THE HARTREE-FOCK ANALYZE 

OF THE STABILITY OF THE 
SUPERCONDUCTING, ANTIFERROMAGNETIC 
AND FERROMAGNETIC STATE 

The free energy at T = of the superconducting state 
for the model (|jj) is expressed by 

n 2 

F s /N = -pt HF + 8XA Ax + UA 2 + U— , (CI) 
where 

t HF = 2t-2Xn, p = 5^(4^) 

A =< CiiCi T > , Ai =< c a c iT > . (C2) 
The stability conditions for this state are: 
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N ^ K E k > 1 

k K 



An 



AT ^ 



2^ fe ' 

Vi + V2Vk 
2E k 



(C3) 



(C4) 



(C5) 



(C6) 



We use the following notation: e k = —tHFVk — Vi = 
4XAi + UA , V 2 = 4XA 0) E k = y/e* + (Vi + V 2Vk ) 2 
and [i is the chemical potential. In general, the solution 
of these self-consistent equations is performed numeri- 
cally. The critical value Us for the stable superconduct- 
ing solution may be found expanding the integrals in a 
series of small parameters An and Ai. The results for 
the rectangular DOS is 

W = 2*\l-Xn/t\ ^ " ^ + 2X + ^ ~ ^ • 

(C7) 

(see also Ref.[ go| for the other DOS). 

The free energy of the AF state in our model is 

o 

F AF /N = -pt HF - Urn 2 + U— , 



(C8) 



where m is the local magnetic moment. It is similar to 
that one in the Hubbard model with the renormalized 
hopping integral tn F . For n = 1 and the rectangular 
DOS it may be expressed in the analytical form as 

with the magnetic moment 

W\l-X/t\ 



(CIO) 



2t/sinh( w ' 1 - x / t ' ) ■ 
The free energy of the ferromagnetic ordering is 

F F /N = -pt HF - AXm{p ] - pi)- Um 2 + ^-n 2 , 

(Cll) 

where Per = {c\ a Cj a ). For the rectangular DOS the stable 
solution is for U > U F , where 

U F _ |1 - Xn/t\(2 - n) - 2(1 - n) 
W ~ n ' 

The value of the magnetic moment is then maximal m = 
n/2 (the saturated ferromagnetism) and the free energy 
is 



(C12) 



W 

F F /N = -—n{\ 



(C13) 



We assumed n < 1. The results f or n > 1 one obtains by 
the proper rescaling. The value (C13) is, of course, the 
same as the free energy for the singly occupied normal 
state obtained in the SBMFA. 



APPENDIX D: THE ANTIFERROMAGETISM IN 
THE HUBBARD MODEL WITH CORRELATED 
HOPPING ANALYZED BY MEANS OF THE 
SBMFA 

We want to find the free energy of the antiferromag- 
netic state in our model, and in order to do it we need 
first to formulate the Hamiltonian in the SBMFA for this 
state. The spin operator Si is given by 



•Pi 



(Dl) 



On a hypercubic lattice the magnetic moment is modu- 
lated with the wavevector Q = (7r, 7r, ...) and its absolute 
value is given by 



in 



= E <J (4a^) exp(iQ • R*) 
= E a (pLoPi™) exp(iQ • Ri) 



(D2) 



where Ri is the lattice vector. To express the hopping 
term of the Hamiltonian we use the renormalized Hub- 
bard operators in the form (^|) neglecting off-diagonal 
elements of the operators p, p, b and b. The fermionic 
part of the Hamiltonian is 

H = Y^K-MAFVk + A ) ]T fljka 
fc cr 

+A m ^a/t CT / fc+Qff ] , (D3) 

where the effective hopping integral 

b 2 [VI - b 2 + 2m + \/l-b 2 - 2m} 2 



t AF = {t-X) 



The free energy is given by 



1 - 4m 2 



F/N = - 1 J2 E k + U b - - 2A m m 



(D4) 



(D5) 



where E k — ^ At 2 AF r] 2 -\- \ 2 n . We used the fact that 
An = for n = 1. The stable solution is found from 
the minimum of F, i.e. from the conditions dF/db = 0, 
dF/dm = and dF/d\ m = 0. We see that the corre- 
lated hopping changes only the effective hopping integral 
tAF-, and therefore, the result for any X is that one for 
X = with the proper rescaling. 
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FIG. 1. Free energy difference of the normal and the 
superconducting state calculated with in the SBMFA at 
n = 1, 0.8, 0.5, 0.2 (from top) for the ID Hubbard model with 
correlated hopping at X = t. It is shown the value of b 2 in 
the normal phase. 



FIG. 2. Free energy at T = for the ID Hubbard model 
with correlated hopping at X = t, n = 1, calculated by means 
of the SBMFA (solid curve)nthe HFA (long-short dashed 
curve) and the exact method [u] (dashed curve). Insert shows 



the differences F ex — F sb (solid curve) and F e . 
curve). (The bandwidth W = 4t.) 



-F Ht (dashed 



FIG. 3. Stability diagram for the superconducting state of 
the Hubbard model with correlated hopping at X = t deter- 
mined by the SBMFA for the rectangular DOS (short dashed 
curve) and the ID chain (solid curve), the exact method for 
the ID chain (curve with dots), and by the HFA for the rect- 
angular DOS (long dashed curve) and for the ID chain (curve 
with crosses). 



FIG. 4. Stability diagram for the superconducting phase 
of the Hubbard model with correlated hopping at n — 1 de- 
termined for the rectangular DOS by the SBMFA (below the 
solid curve). The solutions for the normal phase are: the Mott 
insulator (above the dashed curve), the diamagnetic metal 
(below the dotted curve) and the correlated metal is between 
them. For comparison, the superconducting solutions in the 
HFA are stable below the long-short dashed curve. 
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FIG. 5. Free energy difference of the normal and supercon- 
ducting state as a function of U for the rectangular DOS at 
n = 1 and X — t (solid curve), 1.2t (short dashed curve), 
1.5t (long short dashed curve), 0.8t (long short short dashed 
curve), 0.5t (long dashed curve) and —O.St (dotted curve) - 
Fig. (a). Energy gap vs. U for the parameter X — 0.5t (long 
dashed curve) and 1.2t (short dashed curve) - Fig.b. It is 
seen the different dependencies of E g close to the transition 
to metallic phase at Us = 0.2931^ for X — 0.5t and to in- 
sulator at Us — 1.217W for X = 1.2t, respectively. Fig.c 
shows the change of the effective bandwidth W e ff with U 
for the stable solutions at X = 0.5t (long dashed curve) and 
X = 1.2t (short dashed curve). 



FIG. 6. Free energy difference of the normal and super- 
conducting state as a function of X determined within the 
SBMFA for the rectangular DOS at n = 1 and U = 0.3W 
(long dashed curve), U = 0W (solid curve), — 0.3W (short 
dashed curve), — 0.5W (long short dashed curve), —IW (long 
short short dashed curve)- Fig. a. Energy gap E g as a function 
of X at n = 1 and U = -0.3W, 0W, 0.3W (from top) - Fig.b. 
In Fig.c the dependence of W e f / vs. X is shown for the stable 
solutions at U/W = —0.3, 0.0 and 0.3 (from bottom). 

FIG. 7. Stability diagram for the superconducting and the 
antiferromagnetic phase of the Hubbard model with corre- 
lated hopping at n = 1 determined for the rectangular DOS 
by the SBMFA (solid curve). The dashed curve separates the 
normal phases: the Mott insulator (above) and the correlated 
metal (below); the solution for the diamagnetic metal is below 
the dotted curve. For comparison, the boundary line between 
the AF and superconducting phase obtained in the HFA is 
shown by the long-short dashed curve. 

FIG. 8. Stability diagram for the metallic solutions: for the 
strongly correlated (SC) state (between the solid curves), the 
diamagnetic metal (below the dotted curve) and the weakly 
correlated (WC) state (between the dotted and solid curves) 
at the electron concentration n = 0, 0.2, 0.5, 0.8, 0.99 (from 
top) - Fig.a, and for n = 2, 1.8, 1.5, 1.2, 1.01 (from bottom) 
- Fig.b. 

FIG. 9. Dependence of the parameter b 2 on X in the nor- 
mal state for U = -0.75W and 0.25W (from top) at the 
electron concentration n = 0.5 - Fig.a and n = 1.5 - Fig.b. 
6 2 = 2|<5| = |1 — n\ means that doubly occupied sites are for- 
bidden, for b 2 — 1 there are only doubly occupied and empty 
sites. In Fig. c and d the dependence of the effective band- 
width W e ff vs. X for U/W = -0.75 and U/W = 0.25 (from 
bottom) at n — 0.5 and n=1.5, respectively. 



FIG. 10. Stability diagram of the superconducting state for 
the Hubbard model with correlated hopping for n = 0, 0.2, 
0.5, 0.8, 0.99 (solid curve from top for X/t < and from 
bottom for X/t > 0) - Fig.a and for n = 2, 1.8, 1.5, 1.2, 
1.01 -(solid curve from top for X < Q.5t and from bottom for 
X > 0.5) - Fig.b. The broken curve is a boundary between 
the normal states of strongly and weakly correlated electrons, 
below the dotted curve there are solutions for the diamagnetic 
metal (as in Fig. 8a and Fig. 8b). 



FIG. 11. Comparison of the stability diagram for the su- 
perconducting state obtained in SBMFA (solid curves as in 
Fig. 10) and HFA (long-short dashed curves) for n = 0, 0.2, 
0.5, 0.8, 0.99 (solid curves - from top for X/t < and from 
bottom for X/t > 0; long-short dashed curves - from top for 
X/t < and from bottom for 1/n > X/t > 0, two curves for 
X/t > 1/n correspond to n = 0.99 and 0.8) - Fig.a and for 
n = 2, 1.8, 1.5, 1.2, 1.01 (from bottom for X/t < 0, from top 
for X/t > and from bottom for X above the singular point) 
- Fig.b. 



FIG. 12. Dependence of the gap E g in excitation spectrum 
on the electron concentration n obtained in the SBMFA (solid 
curves) and the HFA (long short dashed curve) for X/t — 0.4, 
U/W = -0.2 and X/t = 0.4, U/W = 0.2 (curves from top 
and the left hand side) and X/t = -0.8, U/W = -0.2 and 
X/t = -0.8, U/W = 0.2 (curves from top on the left hand 
side) . 



FIG. 13. The gap E g in the excitation spectrum (Fig.a) and 
the free energy difference of the normal and superconducting 
state (Fig.b) as a function of n determined within the SBMFA 
for X = 1.2t and U = 0.2W (solid curve), and for X = 3t 
and U = 4W (dashed curve). The values of the dashed curve 
in Fig.a are 10 times smaller. 
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